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Superstring theory predicts the potential formation of string networks with bound states ending 
in junctions. Kinematic constraints for junction formation have been derived within the Nambu- 
Goto thin string approximation. Here we test these constraints numerically in the framework of 
the Abelian-Higgs model in the Type-I regime and report on good agreement with the analytical 
predictions. We also demonstrate that strings can effectively pass through each other when they 
meet at speeds slightly above the critical velocity permitting bound state formation. This is due to 
reconnection effects that are beyond the scope of the Nambu-Goto approximation. 
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I. INTRODUCTION 

There has been a revival of interest in cosmic strings 
since it was realised that many fundamental string the- 
ory models predict so-called cosmic superstrings. The 
leap from string scale to cosmic dimensions is highly non- 
trivial. At first sight there are many severe problems 
related to the formation and growth of superstrings to 
macrosopic scales, such as too high string tension and 
instability towards breaking up, as pointed out already 
in Subsequent advances in fundamental string the- 
ory have changed this picture. All these issues can be cir- 
cumvented and are naturally evaded in many models that 
possess otherwise desirable aspects from the point of view 
of providing successful phenomenology. For instance, 
in the scenarios based on warped compactifications the 
warp-factor can reduce the string tension. Furthermore, 
in models of brane inflation not only are cosmic strings 
produced at the end of the inflationary epoch 0, 0] , but 
it has been argued in great abundance too [1]. 

Considerations inspired by superstring theory suggest 
a particular kind of cosmic string network that consists of 
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fundamental F-strings, Dirichlet D-strings (or more pre- 
cisely Dl-branes) and bound states of these two, known 
as (p, 5)-strings as an abbreviation for p F-strings and 
q D-strings [5| . The presence of {p, (7)-strings brings an 
additional feature to these networks compared to those 
made of only one type of solitonic strings. Namely, where 
two different types of string meet at a point and form 
a bound state leading away from that point, there is a 
junction in the network. Another property of cosmic su- 
perstrings is that their reconnection probability P can be 
small, P <C 1 (for a review see Q and references therein). 

A substantial body of work has dealt with modelling 
of cosmic string networks, which is a challenging task 
due to the combination of several scales involved and the 
non-linear nature of the problem. For instance, prop- 
erties of networks at small scales are still not entirely 
understood and have been under intense analytic study 
recently 0, H 0]- However, it is well-established that 
string networks lose energy efficiently. A scaling string 
network contributes only a fixed, tiny fraction of the to- 
tal energy budget of the Universe, a property that has 
made cosmic strings viable, in contrast to some other de- 
fect models, and generated the appeal of cosmic strings 
from the very start [lo| . 

For the above mentioned reason the effect of junctions 
and bound states on the evolution of the networks is of 
profound interest: instead of reaching a scaling regime 
the network could end in a frozen-out state and start to 
dominate the energy density of the Universe, which can- 
not be tolerated in the standard cosmologies. To date. 
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the evolution of networks with junctions has already been 
studied in several models. Good evidence for scaling was 
reported in [llj with an SU(2)/Z3 model of global strings, 
the junctions being global monopoles. Another field the- 
ory study used a model involving two sets of U(l) 
gauge fields (see also [l3|). There bound states are re- 
ported to have a significant effect on the network in the 
absence of long-range interactions, whereas in the case of 
global strings j unctions play a minor role. Other stud- 
ies include [3,13, EB], which concluded that, even being 
conservative, the networks can scale also in the presence 
of junctions. 

If the bound states and junctions are of major impor- 
tance, then a natural question to pose is under which con- 
ditions they form. This was examined analytically based 
on the Nambu-Goto action in [13, It was shown that 
strong kinematic constraints apply to the formation of 
the bound state. The purpose of this study is to test these 
constraints in a field theory set-up. We work within the 
Abelian-Higgs model in the type-I regime, investigating 
when the intersection of strings results in the formation 
of a bound state (also called a zipper). This topic has 
been addressed already before both analytically [lOj and 
with numerical experiments [20| . Here we revisit it due to 
the renewed interest in more complicated networks with a 
spectrum of string tensions and with improved computa- 
tional resources available. Before presenting the results, 
we briefly review the outcome from the Nambu-Goto ap- 
proximation and introduce the model together with the 
numerical approach. 



II. STRING JUNCTIONS 

Consider a straight string making an angle a with the 
positive X-axis and another one with an angle —a (the 
total angle between the strings being thus 2a). Both 
strings are on the xy-plane and have a velocity v along the 
z-axis with opposite directions and have string tensions 
fii and /i2, respectively. Once they intersect, these can 
potentially form a third string, an "a;-link" which has 
tension (we follow here the notation introduced in ; 
if fii — fi2 then the x-link is indeed positioned along the 
X-axis) . 

The approach in [IE] is based on studying the ac- 
tion at a junction where three strings meet. The kine- 
matic constraints follow from the requirement that the 
total length of the progeny string must increase. The 
allowed parameter region for the link formation can be 
determined at least numerically for any combination of 
string tensions [l^. However, when the colliding strings 
have the same tension, denoted hereafter by fii , it is pos- 
sible to express this requirement in a simple closed form 

a < arccos (- — ) , (1) 



where 7 = — v'^ and necessarily fig, < 2fii. 



Here we want to test this result numerically in a field 
theory where strings are solitonic objects with a non- 
zero width and see if bound state formation is the dy- 
namically preferred process. 

III. MODEL AND NUMERICAL 
IMPLEMENTATION 

The Abelian-Higgs model is governed by the La- 
grangean 

The model has vortex solutions [2lj, in which the scalar 
field can be expressed, with the help of a function / of 
the radial distance r only, as follows: 

= /(Oe™', 

where n is the winding number. For n — 1 the vor- 
tex solutions are topologically stable for any value of the 
scaled couphng, (3 = \/2q^ = "^scaiar/'^gaugo where q is 
the gauge coupling, mgcaiar = VA?? and Wgaugo = \/2gri. 
The domain /3 < 1 is known as the type-I regime. There 
a string with winding number n has lower energy than 
n strings with winding number n = 1 together. In par- 
ticular, two n = 1 strings can merge to form an n = 2 
string [22]. In other words, using the previous notation, 
Us < 2/11 , and thus this regime provides a testing ground 
for the formation of links. 

The numerical code used to evolve two boosted strings 
on a lattice has been reported in [2^ . The boundary con- 
dition employed here was introduced in [13] : the fields on 
the boundaries are updated as though the strings moved 
undisturbed at constant, initial velocities. Of course, this 
means that the simulations cannot be trusted much be- 
yond the time when the kinks on the strings generated 
by their interaction reach the boundary. The majority 
of the simulations were carried out in a computational 
grid of size 400"^, setting the lattice spacing dx to be 0.2 
(the physical size of the lattice in linear dimension being 
therefore 80) and the time step dt = 0.065 {dx and dt 
here in units of {qr])~^). The robustness of the results 
has been tested by varying the size of the simulation box, 
lattice spacing and the initial separation of strings. All 
the snapshots presented in the following section are from 
simulations performed on a 600'^ grid, setting dx = 0.15. 

IV. RESULTS 

Local cosmic strings with an equal winding number, 
here n = 1, always intercommute when intersecting. 
However, after intercommutation, strings can come to- 
gether and merge to form a string of winding number 
n = 2. This occurs when strings meet at relatively small 
velocity and angle (see [l^, [lo]). Alternatively strings 



FIG. 1: Snapshots (from left to right and top to bottom) showing constant energy density isosurfaces in a simulation of two 
strings with n = 1 at /3 = 0.125 colliding with velocity v — 0.1 (top one moving downwards and vice versa), at an angle 
a = 15 °, and forming an x-link as indicated by the o in Fig. (2] 




FIG. 2: Kinematic constraints in the (a, ?;)-plane for jS = 0.36 (left) and l3 — 0.125 (right). The solid line shows the curve of 
equality in (IT|, + marks the events of single intercommutation, whereas o indicates strings merging together to form an x-link, 
a string with winding number n — 2. Events where strings pass through each other by intercommuting twice are shown by A. 
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FIG. 3; Snapshots (from left to right and top to bottom) showing constant energy density isosurfaces in a simulation of 
two strings with winding n = 1 at /3 = 0.125 colliding with velocity v — 0.5, at an angle a — 10°, and undergoing double 
intercommutation, as indicated by the A in right hand Fig. [2] 




FIG. 4: Snapshots (from left to right) showing constant energy density isosurfaces in a simulation of strings with winding n — 1 
(initially above) and n = 2 at f3 — 0.125 colliding with velocity v = 0.1 and at an angle q = 10 °. 
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only intercommute and never meet again. In Figure [T] 
four snapshots are presented that show the isosurfaces of 
constant energy density (set to be 0.2 in dimensionless 
units, roughly 30% of the maximum value in the core for 
a string with n — 1, a, value chosen so that strings with 
winding rt = 1, 2 and 3 can all be readily identified) when 
two strings with n — I collide at velocity v — 0.1 and an- 
gle a = 15° and the coupling constant is /3 = 0.125. 

The formation of the bound state in the Abelian-Higgs 
model is not entirely a straightforward, instantaneous 
process, as already reported in [23|. The link forms 
and unforms, the two strings becoming separately visi- 
ble again. This is not surprising; the process is essen- 
tially right-angle scattering which is well-established by 
numerical experiments with vortices [25l [2^ . Analyti- 
cally, right-angle scattering after a head-on collision is 
understood at critical couphng {(3 = 1) 27], in terms of 
geodesic motion in the moduli space approximation. This 
approach was introduced for the first time in the con- 
text of monopole scattering [2^ but has been recently 
re-employcd in studies of strings in more complicated 
models [29, 30,]. After intercommutation, here the string 
segments align almost parallel before coming together. 
Therefore locally the second collision takes place head-on 
with effectively zero impact parameter and strings scatter 
perpendicularly with respect to the incoming direction. 
Repeated right-angle scatterings are clearly visible in the 
snapshots; the time-scale of attenuation is comparable to 
the time the simulations can be evolved. 

The main result of this study is presented in Figure [21 
where the results from the simulations, standard inter- 
commutation versus link-formation, are shown together 
with the prediction of ([1]) for the allowed parameter re- 
gion in the (a, u) -plane. A decrease in the scalar coupling 
A reduces the tension of a string with winding number 
n = 2 relatively more than that of a string with n — 1. 
This allows one to vary the ratio iJ-s/fii to a certain ex- 
tent by reducing the parameter [3] a.t [3 — 0.36 we obtain 
^sZ/ii ~ 1.88, whereas fi = 0.125 yields fJ-s/ni — 1.77. 
This already leads to a considerable difference when the 
limiting curve of the equality in ^ is plotted for /3 = 0.36 
(left) and /3 = 0.125 (right) in Figure [3 Standard in- 
tercommutation events are indicated by crosses and the 
formation of links by circles. At (3 = 0.36 the agreement 
between the simulations and the analytical prediction is 
perfect within the precision of the grid used in the (a, v)- 
parameter space. There is a discrepancy at the largest 
angles when (3 = 0.125; while the formation of an x-link 
should be allowed, this is not observed in the simulations. 

This behaviour is to be expected from considerations 
of energy conservation. In the Nambu-Goto model, all 
the energy released by shortening of the colliding strings 
goes into the formation of the linking string. In the field- 
theory model, by contrast, some energy is radiated away. 
We therefore expect that it should be slightly harder to 
form a link, and so there should be a small band where 
link formation is possible for Nambu-Goto strings, but 
not for field-theory strings. 



In addition, at /? = 0.125 there is a band just above 
the highest velocity allowing x-links to form suggested 
by ^ where a link does not form, but strings come to- 
gether again and intercommute for a second time, events 
denoted by triangles in Figure [2] Snapshots of this pro- 
cess are presented in Figure [H The end result is thus 
two strings consisting of the same segments as initially 
and indistinguishable from the situation where strings 
had passed through each other apart from some defor- 
mation around the interaction point. A similar type of 
effective non-intercommutation has been reported when 
strings collide at very high velocities [23l . [2^ ] . There is no 
sign of this kind of phenomenon at /? = 0.36. This may 
be because less energy is radiated away in the field when 
the scalar mass is larger. 

The general case, where all three strings have different 
tensions can be investigated in the Abelian-Higgs model 
by colliding strings with different winding numbers. This 
was done for two strings with windings n = 2 and n — 
1 at /3 = 0.125. Though the explicit analytic formula 
for the asymmetric case (sij is not presented here, the 
kinematically allowed area is almost degenerate with the 
one presented in Figure [2] for the coupling j3 = 0.125. A 
complication in the Abelian-Higgs model is that due to 
unequal winding numbers, the intercommutation leads to 
a formation of a bridge between the strings (see also [s^). 
whose influence on bound-state formation at the very 
least cannot be entirely neglected. No systematic study 
in the kinematic parameter space was carried out, but 
Figure [4] shows snapshots of a simulation that provides 
evidence for a bound state formation now in the form of 
a string with winding number rt = 3. 



V. DISCUSSION 

We have reported on numerical simulations of string 
collisions in the Abelian-Higgs model. The objective 
was to monitor the kinematic parameter space (a, v) and 
compare the outcome to the analytical prediction. We 
do not observe x-links outside the area where they are 
not expected to be kinematically allowed. On the other 
hand, bound states generically form whenever allowed 
to appear. This is interesting because such links do not 
have to form dynamically; there could have been simple 
intercommutation events instead, but almost the whole 
region appears to prefer to form x-links. 

It is not surprising that the observed discrepancy oc- 
curs when strings intersect at large angles - the Nambu- 
Goto action does not include the effects of intercommu- 
tation. As strings in the Abelian-Higgs model always in- 
tercommute, effectively the original strings break up, and 
when the new strings straighten after this reconnection, 
they do not come together and merge at large angles. 
This could be different for strings that do not intercom- 
mute and exchange partners and it would be interesting 
to see what a study with strings in two separate gauge 
fields would yield. 
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Secondly, we have demonstrated once more that the 
'effective' reconnection probabUity of strings even in the 
Abehan-Higgs model is not strictly unity. This is well- 
documented by numerical experiments at high collision 
speeds where a second intercommutation even takes place 
(see e.g. [l^ but as shown here can occur at mod- 
erate velocities too (it would be interesting to see if a 
modification of the moduli space approximation could 
capture the second intercommutation event). It was ar- 
gued in (20| that when the bound state with winding 
n = 2 dissolves (which inevitably happens eventually 
with the boundary conditions introduced in [IJ), the 
original strings re-emerge. This is effectively like no in- 
tercommutation taking place. However, the process re- 
ported here does not proceed via bound state formation 
and dissolution, but rather by repeated intercommuta- 
tion, and seems therefore to be of a different nature. 
This is further confirmed by the inequality ([1]) accord- 
ing to which the bound state is forbidden at velocities 
where the double intercommutation is observed. 

Obviously the use of the Abelian-Higgs model has re- 
stricted us in this study to very limited values of the 
ratio of string tensions fi^/^i, whereas models inspired 
by superstring theory would allow a much wider range of 
values. We have shown that strings in the Abelian-Higgs 
model can have the same principal features as expected 



from cosmic superstrings, namely formation of junctions 
and effective reconnection probability not strictly unity 
even at moderately low velocities. While network simu- 
lations very deep in the Type-I regime are hardly feasible 
to perform, these results may have implications for theo- 
ries relying on a very low value of the scalar coupling A. 
Such a network can arise e.g. upon breaking a gauge sym- 
metry along a fiat direction in supersymmetric theories 
as has been pointed out recently in [33|. To summarize, 
as anticipated already in [T^ . this study demonstrates 
the richness of strings even in the simple Abelian-Higgs 
model. 
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